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Abstract 

We study the 2-flavor lattice Schwinger model with Wilson fermions 
in the chiral limit. The quark mass is determined using the PCAC 
definition. We numerically compute the masses of the iso-triplet (vr) 
and iso-singlet particles (ij) for diff'erent quark masses and compare 
our results with analytical formulas. 
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Introductory remarks 

Although we begin to understand the formulation of chiral fermions on the 
lattice [I], the high cost of their simulation still renders Wilson fermions 
an attractive alternative. Wilson fermions, however, suffer from an explicit 
breaking of chiral symmetry and can only in a carefully taken chiral limit 
be used to analyze chiral aspects of QCD. Approaching and properly under- 
standing this limit is a non-trivial task and studies in 2-dimensional models 
contribute to our understanding. 

The Schwinger model j2] with two flavors is a particularly convenient test- 
bed for studying the chiral limit since it shares several characteristic features 
with QCD. The low lying spectrum ("mesons") contains an iso-triplet of 
light particles (tt's) which become massless in the chiral limit. In addition a 
massive iso-singlet state {rj) appears, giving rise to a U(l) problem equiva- 
lent to QCD. For the massless case the analytical solution of the Schwinger 
model in the continuum has been known for many years and also for small 
fermion masses some analytical results for the model with flavor have become 
available ^3j-^12.. 

For the lattice model with staggered fermions the connection to these 
continuum results is well understood |ISI-|2D|- The challenge in a lattice 
simulation with Wilson fermions is to properly fine-tune the bare parameters 
and drive the pions to zero mass. Although many aspects of the Schwinger 
model with Wilson fermions are well understood - in particular the relation 
between the spectrum of the Dirac operator and the topology of gauge fields 
has been studied intensively - relatively little I^HJ-EHl, [U] is known 
for the chiral limit with Wilson fermions. What is still lacking is a high 
precision numerical study of the mass spectrum for Wilson fermions com- 
bined with a systematic comparison to analytical results. Furthermore the 
Schwinger model allows for a simple implementation of the PCAC method 
pUj for defining the quark mass. Using this definition together with the nu- 
merical results for the mesons in the comparison to analytical formulas leads 
to an interesting test of the PCAC quark mass definition. 

In this letter we report on our study of the two-flavor lattice Schwinger 
model with Wilson fermions. In particular we concentrate on the low-lying 
particle spectrum (vr's, rj) and the behavior in the chiral limit. We use the 
PCAC method to define the quark mass m and compare our results for the 
spectrum of the meson masses as a function of m with analytical formulas. 
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Setting, quark mass and observables 



The action for the Schwinger model with Wilson fermions consists of two 
parts S = Sg + Sp with Sq and Sp being the gauge field and fermion action 
(2 flavors) both in standard Wilson formulation. We work on L x L lattices 
with periodic boundary conditions for the gauge fields and mixed periodic 
boundary conditions (periodic in space-, anti-periodic in time-direction) for 
the fermions. We simulate the system on lattices of size L = 16, L = 24 
and L = 32 and /3-values of /3 = 2, 4, 6. For the smallest lattice we use 10^ 
configurations for the determination of the observables, while for L = 24 and 
L = 32 we typically use 10^ configurations but increase this number for some 
values of the parameters where better statistics is desirable. The update 
is done using the standard hybrid Monte Carlo method |22| with 10-steps 
trajectories and a step size adjusted such that the average acceptance rate 
is 0.8. Several observables, including the geometric topological charge, are 
monitored and appropriate numbers of configurations are discarded between 
the measurements to assure proper de-correlation and achieve topological 
ergodicity. 

The quark mass is determined using the PCAC method as outlined in pUj . 
The essential idea is to use the coefficient of the right hand side of the axial 
SU(2) Ward identity as a definition of the mass. One computes the matrix 
elements with a pseudo-scalar source on both sides of the Ward identity and 
the quark mass is then given by the ratio of the vev.'s for the two sides. The 
details on the implementation of the method for the lattice Schwinger model 
were reported elsewhere [21] • Here we just briefiy summarize the relevant 
formulas. The quark mass is defined as (x 7^ /i = 1, 2) 



m 



1 f(o|<t /^lo) - {oiK^tW 



(0|%Vs|0) 



(1) 



with 



i'x ® 75 , 



The gamma matrices are chosen as 71 = (Xi , 72 = o"2 and 75 = with aj 
denoting the Pauli matrices. The generators of rotations in fiavor space 
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are also given by the Pauli matrices r" = (Tq, a = 1, 2, 3. The quark mass m 
was computed for each /3, L and k individually. 

The masses for the pions and the eta-particle were computed from the 
decay of two point functions of the following currents 

JT = V'x ^'^ ® 7m for the vr's , 

= l®7^Vx for the r/. (3) 

We emphasize that 2-point functions of scalars or pseudo-scalars are not 
particularly well suited for the determination of the meson masses. These 
operators are bosonized (compare the discussion below) by cosines and sines 
of the fundamental fields [HI EI and their 2-point functions strongly mix con- 
tributions from both the triplet and singlet states. 

Finally let us briefly discuss the correspondence of the gauge coupling on 
the lattice to the gauge coupling g as it appears in the continuum formulas 
below. The Schwinger model is a super-renormalizable theory, i.e. the bare 
coupling does not get renormalized and equals the physical coupling g. 

On the lattice we use the usual factor P in front of the gauge field action. 
In the naive continuum limit it is related to the continuum coupling g via 
(we set the lattice spacing a = 1) 

Again we can use the super-renormalizability of the model and define the 
coupling g which we use for comparison with the analytical results in the 
continuum (see below) through (jlj. 



Analytical results for the continuum model 

Before we start with the presentation of our numerical results, let us briefly 
discuss the known analytical results for the continuum model. This is best 
done by looking at the generalized Sine Gordon model which bosonizes the 
2- flavor Schwinger model [3]. Its Lagrangian reads 

l(5y,{0))^ + bdy,^^)? + i ^U^)? - cos (v^ cos (v^¥,(3)) . 

£j ^ TT 7T 

(5) 

The constant c is not determined by the bosonization but is related (see 
e.g. |S1 CHj) to the masses /i^^-* and [i^^^ used for normal ordering of the 
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fields v?'-"^ and (^9*^^) by c = (fi^'^^ fi^'^^Y^'^e^ /2, where 7 denotes Euler's constant 
7 = 0.577216 .... For the case of two flavors two of the currents which we 
measure (compare Q) are bosonized with the following prescription (here 
rO = 1) 

?^(x) ® 7^ ^'(a;) = -^e^,d,if^''\x) , a = 0,3. (6) 



For the other two members of the iso-triplet no abelian bosonization 

is known. However since the original model is invariant under flavor rotations, 
their masses are the same as for the triplet current (a = 3, pions). When 
setting the quark mass m in (jSJ to zero, the two flavor Schwinger model is 
bosonized by two free fields. One of them, the pion field is massless 
(Mtt = 0) and the eta field <y9*^°^ obtains the Schwinger mass = g^2/TT. 

For non-vanishing quark mass, also the bosonized model © can no longer 
be solved in closed form. A first approach is to use a semi-classical analysis 
(see e.g. ^U]) which gives a good approximation when all involved masses 
are large. One simply determines the minimum of the interaction part 
of © and the squares of the masses are then given by the second 
derivatives of ^p^^^] at the minimum. After normal ordering the fields 

with respect to their own masses (i.e. setting jj^^^ = and n^^^ = M^^), 
the semi-classical analysis for the iso-triplet and iso-singlet gives 



g IT' 



and 




9 V ^ g ' 
We remark that the latter equation, i.e. the semi-classical relation between 
and M^, holds independently of any normal ordering prescription or the 
value of M^, due to the formula d'^/d^p^'^^ V = g'^2 / ir + d'^ / d(p^^^ V, which in 
turn is evident from (jS)). This will be explored below. 

In an attempt to go beyond semi-classical analysis one can approximate 
the generalized Sine-Gordon model © by a solvable model. One possibility 
|1] is to analyze the model in the limit of large coupling g and small mass m 
where the iso-singlet field (p^^^ becomes static and the model is reduced to a 
standard Sine-Gordon model for the triplet field (f^^^ 



(3)r _ 2Ccos(v^v5(='M . (9) 
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This reduced bosonized theory has been studied jH HHI using the WKB ap- 
proximation [nni for (jni). The spectrum of the model (jHI) contains a sohton, 
an anti-sohton and two sohton-anti-sohton bound-states. The hghter bound 
state has a mass equal to the soliton (and anti-soliton), and the mass for this 
triplet can be related to the mass of the iso-triplet in the 2-flavor Schwinger 
model. More recently © was studied again by Smilga [12] but his analy- 
sis is now based on the newly found analytic solution [HI] of the standard 
Sine-Gordon model. Smilga finds 





2.008... — . (10) 



Obviously the truncation of the original model © to the standard form 
does not allow for a result for the singlet mass M^. For this state only the 
semi-classical formula (jH)) is available. It his however interesting to use the 
exact result (fTUj) of the truncated model as an input in (jH)), and below we 
compare also this formula to the numerical data. 



Numerical results 

Let us now discuss our numerical results for the mass spectrum. We start 
with presenting the data for the triplet mass in Fig. Q and compare them to 
the semi-classical-result (0) and Smilga's formula (fTIH) . We show our results 
for /5 = 2, 4, 6 in the first three plots and and in each plot combine all available 
values for the sizes L. In the last plot we zoom in on the chiral region of 
the /? = 4 plot. We compare the numerical results to Smilga's formula (jlOp 
(full line) and the semi-classical result ((7j) (dotted line). We remark, that the 
errors in the determination of m from (P) are smaller than the symbols and 
thus no horizontal error bars were included in the plots. 

It is obvious, that the approach to the chiral limit is under good control, 
i.e. the pion mass M^^ vanishes as the PCAC quark mass goes to zero. In 
addition the data is well described by the analytical result (fTUj) (full line). 
In particular at small quark mass and low (3 where Smilga's assumption 
m -C is fulfilled best, the coincidence of the data with pO|l is convincing. 
For increasing quark mass we observe a cross-over behavior and the data is 
then better described by the semi-classical result ((7j) which is represented by 
the dashed line. This deviation from ()10|) was also observed for staggered 
fermions in [T^. The comparison of the data at different lattice sizes L 
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shows, that the data points essentially fall on top of each other and for our 
parameters L and (3 the finite size effects are negligible. When zooming in 
on the very small m region (last plot in Fig. the error bars become larger 
due to the small M^r, but the approach to the chiral limit is still relatively 
clean. 
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Figure 1: Results for the pion mass M^^. Symbols: Monte Carlo data. Full 
line: Smilga's formula (fTUI) . Dashed line: the semi-classical formula ((Tj). 

In Fig. 121 we show our results for the eta mass at /3 = 4 and 6. We 
find that the data approaches the correct m = value (dotted line) as the 
quark mass vanishes. The approach is smoother as one goes closer to the 
continuum limit, i.e. increases (3. The semi-classical formula (jH)) provides a 



6 



reasonable description of the data for larger m, in particular when Smilga's 
result (fTn|l for M^, is used as an input in (jHl) (full curve). 

For small m the quantum fluctuations become larger and thus the nu- 
merical data deviates from the semi-classical curve. This behavior at small 
m can be fitted by a power law as has been done in ^H] for the results from 
the model with staggered fermions. 
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Figure 2: Results for the eta mass M^. Symbols: Monte Carlo data. Dashed 
curve: the semi- classical formula with taken from (|7j). Full curve: the 
semi-classical formula (jHl) with M^^ taken from (fTIUl . The dotted horizontal 
line corresponds to at quark mass m = 0. 



To summarize, we find that the approach to the chiral limit is under good 
control. As the PCAC quark mass m goes to zero we find that the pion mass 
vanishes and approaches the Schwinger value g^2/'K. The data for 

is well described by Smilga's formula (fTIH) for m g with a cross-over to 
the semi-classical formula (|7j) as m is increased. For the eta mass the data 
is reasonably well described by the semi-classical formula (jHl) with quantum 
fluctuations becoming less important as the quarks are made heavier. 

Acknowledgement : We want to thank Andrei Smilga and Renate Tepp- 
ner for discussions. 



7 



References 



[1] F. Niedermayer, Nucl. Phys. B (Proc. Suppl.) 73 (1999) 105; 
H. Neuberger, Chiral fermions on the lattice, Plenary 
talk at Lattice-99, Pisa, Italy, June 29tli - July 3rd 1999, 



liep-lat/99090042; 



M. Liischer, Chiral gauge theories on the lattice with exact gauge 
invariance, Plenary talk at Lattice-99, Pisa, Italy, June 29th - July 
3rd 1999. 

J. Schwinger, Phys. Rev. 128 (1962) 2425. 

S. Coleman, R. Jackiw and L. Susskind, Ann. Phys. 93 (1975) 
267. 

S. Coleman, Ann. Phys. 101 (1976) 239. 

L. V. Belvedere, J. A. Swieca, K. D. Rothe and B. Schroer, Nucl. 
Phys. B153 (1979) 112. 

H. Joos, Helv. Phys. Acta 63 (1990) 670; S. I. Azakov and H. Joos, 
Helv. Phys. Acta 67 (1994) 723; I. Sachs and A. Wipf, Helv. Phys. 
Acta 95 (1992) 652. 

M. P. Fry, Phys. Rev. D47 (1993) 2629. 

C. Gattringer and E. Seller, Ann. Phys. 233 (1994) 97. 

J. E. Hetrick, Y. Hosotani and S. Iso, Phys. Lett. B 350 (1995) 
92. 

C. Gattringer, Ann. Phys. 250 (1996) 389 and QED2 
and U(l)-probleni, Thesis, University of Graz 1995 (E-Print 
.hep-th/9503137 ). 

C. Adam, Ann. Phys. 259 (1997) 1. 

A. V. Smilga, Phys. Rev. D55 (1997) 443. 

E. Marinari, G. Parisi and C. Rebbi, Nucl. Phys. B190 (1981) 
734. 



8 



S. R. Carson and R. D. Kenway, Ann. Phys. 166 (1986) 364, Phys. 
Rev. D35 (1987) 1961. 

M. Grady, Phys. Rev. D35 (1987) 1961, Nucl. Phys. B365 (1991) 
699. 

V. Azcoiti, G. Di Carlo, A. Galante, A. F. Grillo and V. Lahena, 
Phys. Rev. D50 (1994) 6994. 

H. Dilger, Phys. Lett. B 294 (1992) 263, Nud. Phys. B434 (1995) 
321. 

C. Gutsfeld, H. A. Kastrup, K. Stergios and J. Westphalen, Nud. 
Phys. Proc. SuppL 63 (1998) 266; C. Gutsfeld, H. A. Kastrup, 
K. Stergios, Mass spectrum and elastic scattering in the massive 
SU{2)f Schwinger model on the lattice, E-print: hep-lat/9904015| 

Ph. de Forcrand, J. E. Hetrick, T. Takaishi and A. J. van der Sijs, 
Nud. Phys. Proc. Suppl. 63 (1998) 679. 

F. Berruto, G. Grignani, G.W. Semenoff and P. Sodano, Phys. 
Rev. D59 (1999) 034504, Ann. Phys. 275 (1999) 254. 

J. C. Vink, Nud. Phys. B307 (1988) 549. 

C. R. Gattringer, I. Hip and C. B. Lang, Nucl. Phys. B508 (1997) 
329, Nud. Phys. Proc. Suppl. 63 (1998) 498. 

P. Hernandez, Nud. Phys. B536 (1998) 345. 

W. Bardeen, A. Duncan, E. Eichten and H. Thacker, Phys. Rev. 
D57 (1998) 3890, Nud. Phys. Proc. Suppl. 63 (1998) 811. 

Ch. Holbling, C. B. Lang and R. Teppner, Nucl. Phys. 
Proc. Suppl. 73 (1999) 936. 

A. C. Irving and J. C. Sexton, Nucl. Phys. Proc. Suppl. 47 (1996) 
679. 

V. Azcoiti, G. Di Carlo, A. Galante, A. F. Grillo and V. Laliena, 
Phys. Rev. D53 (1996) 5069. 



9 



[28] R. Kenna, C. Pinto and J. C. Sexton, E-print hep-lat/9812004[ 



[29] S. Elser and B. Bunk, Nucl. Pliys. Proc. Suppl. 53 (1997) 953. 

[30] M. Bochicchio, L. Maiani, G. Martinelli, G. Rossi and M. Testa, 
Nucl. Phys. B262 (1985) 331; K. Jansen, C. Liu, M. Liischer, H. 
Simma, S. Sint, R. Sommer, P. Weisz and U. Wolff, Phys. Lett. 
B372 (1996) 275. 

[31] I. Hip, C. B. Lang and R. Teppner, NucL Phys. Proc. Suppl. 63 
(1998) 682. 

[32] S. Duane, A. D. Kennedy, B. J. Pendleton and D. Roweth, Phys. 
Lett. 195B (1987) 216. 

[33] R. F. Dashen, B. Hasslacher and A. Neveu, Phys. Rev. Dll (1975) 
3424. 

[34] Al. B. Zamolodchikov, Int. J. Mod. Phys. AlO (1995) 1125. 



10 



